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Abstract 

The characteristics relating to the lateral and longitudinal development of EAS in the energy region of 
10^^ — 10^^ eV have been analyzed in the framework of the QGSJET model and of mass composition of 



> 

^ primary cosmic rays. It is found that at £'o > 5 x 10^^ eV the mean mass composition of primary cosmic 



rays begins to vary as indicated by a rise of (In A) with increasing energy. The maximum value of (In A) 
is observed at £^o ~ (5 — 50) x 10^^ eV. It is confirmed by data of many compact EAS arrays and does 



i-^ not contradict an anomalous diffusion model of cosmic ray propagation in our Galaxy. In the superhigh 

Q energy region (> 10^® eV) the value (In ^4) begins to decrease, i.e. the mass composition becomes lighter 

id 

(/2 and consists of protons and nuclei of He and C. It does not contradict our earlier estimations for the mass 

composition and points to a growing role of the metagalactic component of cosmic rays in the superhigh 

(N 

j> energy region. 
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(S| I. INTRODUCTION gies ~ lO^^ _ iqu ^y ^p ^^ ultra-high energies 

I^ IQis - 1021 eV. 

O *•.-.• J- 1 • 1 • The study of PCR mass composition at En > 

.. An investigation oi superhigh energy cosmic 

•^ rays in the 1950s - 1960s and 1970s - 1990s ^^'' ^^ '' P"^^^^^^ ^^^^ ^^^ *^ ^^^^ '^^'''^'- 

\^ , ,,,,,,• r • 1 •.• • .1 ing extensive air showers (EAS) and by indirect 

3 has led to the discovery oi irregularities in the v / ^ 

1 o 1r^T^ TT way only using the EAS characteristics maxi- 

energy spectrum at the energy ~3xl0eV 

1 i^iq T7 rm [TiTi rr^i 11 r ■ mum scnsitivc to the mass composition. For 

and ~ 10^^ eV yj |2||. The problem of on- ^ 



these purposes, at the Yakutsk complex array 
the characteristics of both longitudinal and ra- 
dial development of EAS are taken. The longi- 
tudinal development of showers is reconstructed 
* s.p.knurenko@ikfia.ysnI5| measuring the EAS Cherenkov light, and the ra- 

dial development measuring the structural func- 



gin of breaks in the energy spectrum requires 
reliable knowledge about the mass composition 
of primary cosmic rays (PCR) from high ener- 
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.ye.sleptsov.iktta.ysn.ru tions of electron and muons with Ethr > 1 GeV. 



In this paper PCR mass composition data ob- 
tained at the Yakutsk EAS array in recent year 
are summarized. 

n. METHOD 

The estimation of PCR mass composition 
was determined by the formula: 
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• Yakutsk, N/F; 

• X^^^, big array; 

• E_ / E (small array) 



■ X^^^, small array 
--Lagutin{2001) 
* Cor:N^-N^ 
o E _ / E (big array) 
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{In A) =J2ai-\nAi 



(1) 



where Oj is the relative portion of nuclei to the 
mass number Ai. In each case the experimen- 
tal data were compared with calculations by 
the QGSJET model for the primary proton and 
iron nucleus in the framework of a superposition 
model: 

pexp _ pp 



(In A) 



X InApe, 



(2) 



pFe _ pp 

where P is a parameter characterizing the lon- 
gitudinal or radial development of EAS, for ex- 
ample, Xmax — a maximum depth of a shower 
development, p^/ps — a portion of muons at the 



fixed distance from a core, N., 



K 



a cor- 



relation between the total number of muon and 
electrons at sea level etc. The energy region of 
10^^ — 10^^ eV, where there are "knee" and "an- 
kle" irregularities in the CR energy spectrum, 
has been considered specially. 

III. RESULTS AND DISCUSSION 

Fig. [1] and Fig. |2] present results on the PCR 
mass composition of the Yakutsk EAS array. 



Figure 1 : Mass composition of CR obtained by dif- 
ferent characteristics of EAS in the region of ultra- 
high energies. 

Data were obtained in the framework of the 
QGSJET model and two-component mass com- 
position (proton-iron nucleus). A few character- 
istics corresponding to the radial and longitudi- 
nal development of EAS [[31 HI |5l [6l| are used in 
this analysis. In each case, the value (In A) was 
determined by an interpolation method accord- 
ing to Q. From Fig.[T] Fig.git follows that the 
mass composition changes reaching the heavier 
in the energy region of(2 — 5)x 10^^ eV and 
then, beginning from 3 x 10^^ eV it becomes 
lighter. 

The curves are calculations according to 
the anomalous diffusion model of cosmic ray 
propagation in the Galaxy (dashed line in 
Fig. [1]) [121 and the two- source model (solid line 
in Fig. |2]) jSl. The calculations have been car- 
ried out separately for each of following groups 
of nuclei: p. He, CNO, N-Si, Fe. The summary 
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Figure 2: Mass composition of cosmic rays at super- 
high and ultrahigh energies. The curve is a calcula- 
tion by Lagutin et al (2001) and (2007) according to 
the anomalous diffusion model for cosmic ray pro- 
pogation. 



It should be noted that in the region after a 
"knee" the estimations of the CR mass compo- 
sition obtained at the compact arrays agree well 
with each other. The same cannot be said of 
the energy region of ~ 10^^ eV (see Fig. |2|, 
where HiRes data point to the faster enriching 
of the primary radiation by light nuclei and pro- 
tons in comparison with the Yakutsk array data, 
which point to the gradual change of the mass 
composition from heavier to light (proton and 
He nuclei, on the whole) at Eq ~ 10^^ eV. 
In both cases, the data point to the existence 
of tendency to the "protonization" of PCRs at 
En>3x 10^*^ eV. 



values of {In A) are shown in Fig. [T] and Fig.|2] 
and for the first case in the Table HI Data in the 
Table|l]point to the change of mass composition. 
In the energy region of 5 x 10^^ — 5 x 10^*^ eV 
the mass composition is heavier than at Eq ^ 
10^^ eV. It contradicts the experimental data 
(Fig. [2]). The calculations [8] according to a sce- 
nario for the particle generation in the sources 
of two different types, assuming the analogous 
model for the CR propagation, showed the pres- 
ence of a fine structure in the dependence of 
(In A) on energy (see Fig. |2]). There are also 
the sharp peaks in our experimental data. In this 
connection, the investigation of the PCR mass 
composition with using a great numbers of FAS 
characteristics and statistical data as possible is 
of particular interest. 



rV. SUMMARY 

A character of energy-dependence of (In A) 
(see Fig. |2]) by the Yakutsk array data points to 
the change of primary particle mass composi- 
tion at the singular points of CR energy spec- 
trum. The {In A) rises with energy after the 
"knee", riching the maximum value (In A) = 
3.5 in the energy interval of (2 — 5) x 10^^ eV, 
and then it begins to decrease. Such a behavior 
does not contradict a hypothesis for the propa- 
gation of CR by the anomalous diffusion laws in 
fractal interstellar medium (Lagutin et al, 2001) 
and according to which the mass composition of 
primary particles varies similar to experimental 
data. At Eq > 10^® eV the value (In A) de- 
creases gradually and at Eq ~ 10^^ eV the mass 



Table I: Calculation results for (In A). 



A H, % He, % CNO, % Ni-Si, % Fe, % (In A) (In A) 



Eo,ey 



experiment 



1 X 10^5 22 21 



20 



1 X 10 



16 



19 19 20 



1 X 10^' 17 17 20 



1 X 10 



18 



18 16 19 



119 



1 X 10'^ 55 17 



1 X 10^0 68 18 6 



17 19 2.14 1.55 ±0.25 



18 23 2.32 2.15 ±0.24 



19 26 2.43 3.16 ±0.23 



19 27 2.44 2.61 ±0.21 



11 1.18 1.66 ±0.18 



3 0.69 



composition is practically represented by He nu- 
clei and protons. This does not contradict cal- 
culations by Berezinsky et al. (91 for the meta- 
galactic model, in which the "ankle" observed in 
the experiments on registration of ultrahigh en- 
ergy cosmic rays can only be formed by the pro- 
ton component arriving from the Extragalaxy. 
Thereby, details of experimental spectrum form, 
for example, "dip", i.e. the decrease of inten- 
sity at Eq ~ 10^^ eV, are most likely caused by 



the interaction of the extragalactic protons with 
relic photons (p ± 7 ^ p ± e~ ± e+). As a di- 
rect evidence of this hypothesis, an anisotropy 
can be used which is associated with the ori- 
gin and sources of cosmic rays. According 
to (lOl [m m, at Eo > 8 X 10^*^ eV the weak 
correlation in the arrival directions of EAS with 
the Galaxy plane and the strong correlation with 
the Extragalaxy plane are observed. The quasars 
can be possible sources of cosmic rays. 
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